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Abstract. A one-dimensional, collisionless plasma given by the Vlasov-Poisson system is considered and 
the stability properties of periodic steady state solutions known as Bernstein-Greene-Kruskal (BGK) waves 
are investigated. Sufficient conditions are determined under which BGK waves are linearly unstable under 
perturbations that share the same period as the equilibria. It is also shown that such solutions cannot 
support a monotonically decreasing particle distribution function. 

PACS. 52.25.Dg Plasma kinetic equations - 02.30.Jr Partial differential equations - 52.35.-g Waves, 
oscillations, and instabilities in plasmas and intense beams 


1 Introduction 

A plasma is a partially or completely ionized gas, in which 
the motion of unbound ions generates electric and mag¬ 
netic fields that strongly affect individual particle motion. 
When a plasma is of low density or the time scales of inter¬ 
est are sufficiently small, it is deemed to be “collisionless”, 
as collisions between particles become infrequent. Many 
examples of collisionless plasmas occur in nature, includ¬ 
ing the solar wind, galactic nebulae, and comet tails. 

One fundamental, one-dimensional model of collision¬ 
less plasma dynamics is given by a system of partial dif¬ 
ferential equations known as the Vlasov-Poisson system: 

dtf + vd^f - Edyf = 0 
d^E = 1 - y fdv. 

Here / = f{t, x, v) represents the distribution of electrons 
in the plasma, while E = E{t, x) is the self-consistent elec¬ 
tric field generated by ions and electrons. The independent 
variables, t > 0 and x, u G K represent time, position, and 
velocity, respectively. Instead of studying a large collection 
of ionic species interacting with the electrons, the density 
of ions is given by a neutralizing background, normalized 
to 1 in the equations above. 
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Of course, the study of plasmas has been an extremely 
rich area within both the physics and applied mathemat¬ 
ics communities for decades. Early work regarding the 
stability and instability of plasmas includes 0, m, m, 
PP] , , and . More recent studies of one-dimensional 

plasma dynamics include 0,0, US], m, and m while 
for general references regarding the V lasov equation and 
models in kinetic theory, such as (VP I and its electromag¬ 


netic counterpart, the Vlasov-Maxwell system, we mention 

0, m, and [24]. 

In the current study we are interested in the stabil¬ 
ity properties of time-independent solutions of the sys¬ 
tem. More specifically, we wish to study well-known steady 
states called BGK waves. In 1957 Bernstein, Greene, and 
Kruskal [4] showed the existence of an infinite family of 


exact solutions to (VP I that have come to be known as 
BGK waves. Since then, the stability or instability of these 
solutions has been of great interest to the plasma physics 
and applied mathematics communities (0, [in], [n], HZ], 
[IS])- In particular, recent experimental identifications of 
electrostatic solitary waves in space plasmas have provided 
further justification for the study the properties of such 
steady states. A BGK wave is a steady solution of the 
form 


f{x,v) = ti + 
E{x) = tj)' {x) 


( 1 . 1 ) 


where </> is a periodic solution of 


(j)"{x) = l- j + (t){x'^ dv. (1.2) 
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In [15] the stability properties of such solutions to (VP I 
were studied. There it was determined that any BGK wave 
is linearly unstable with respect to multi-periodic pertur¬ 
bations, that is, pertnrbations whose periods are integer 
multiples (2 or greater) of the period of the steady state (f>. 
Unfortunately, such a general result in the singly-periodic 
case, in which the perturbation and BGK wave possess 
the same period, remains unknown and does not seem 
likely. As the author states in [15], one cannot hope that 
all periodic BGK waves are unstable with respect to such 
perturbations. Instead, the stability properties may de¬ 
pend delicately on the specific BGK wave of interest. To 
date little information regarding this question is known, 
though some results concerning weakly-inhomogeneous so¬ 
lutions have been obtained m, m)- The main focus of 
our paper, then, is to determine sufficient conditions un¬ 
der which these BGK waves (/r, (j)) are linearly unstable 
due to singly-periodic perturbations. 

As Lin constructed a strong framework for studying 
the problem, we will follow many of the assumptions and 
previous results developed in m- To be precise, we make 
the following assumptions on the BGK waves dull ) - dL^ 
throughout: 

(i) fJ, S C^(R) is nonnegative. 

(ii) /r satisfies the condition of neutrality, i.e. 

/ = 

(iii) fi' decays at infinity, i.e. there is 7 > 1 and C > 0 
such that 

(iv) We arrange the period of in a specific manner. Let 

be the minimal period of the solution cj) from 0 
so that 

(l){x) = (l){x d- P^) 

for every a; G K and define the quantities 

6- = min d>(x) 
a:G[0.P,^] ' 

6+ = max (bix). 
xe[o,P4 


Then, without loss of generality we may rearrange 
the starting point to impose conditions on the values 
of (j) so that it satisfies: 


</'(0) = (b{P<t,) = (b+, 



cbix) = cj){Ptp - x), Va; G [0, P^]. 


Additionally, without loss of generality we can take 
(j) to be strictly decreasing on the interval 0 , 


With this structure in place, we will show that (0 
cannot possess solutions for which ^ is strictly decreas¬ 
ing. More importantly, we will show in subsequent sec¬ 
tions that, while one cannot hope that all such solutions 


to ( 1 . 1 )-( 1 . 2 ) are linearly unstable with respect to P^- 
periodic perturbations, conditions do exist which guaran¬ 
tee their instability at the linear level. To be more precise 
further definitions are needed. For p, given above we define 
the function 


q{x) = j p! + (bix'^ dv 

and let Aq < 0 be the smallest eigenvalue, with corre¬ 
sponding first eigenfunction ipQ G H^{0, P^,), of the Sturm- 
Liouville [5^ problem 

tb"(x) + {q{x) -f \)tb{x) = 0 

with ib{0) = tbiP,j>) = 0. We will prove that if q satisfies 
J {q{x) + Ao)(g(a:) -k ^Ao)|V'o(a;)|^ dx <0 (v) 

the BGK wave (/i, (b) must be linearly unstable with re¬ 
spect to perturbations of period P^. In particular, this 
demonstrates that for a BGK wave to be stable, the val¬ 
ues of q cannot be essentially localized within the interval 
[—Ao,—|A o]. Our main results will be stated precisely in 
the next section, after deriving the specific partial differ¬ 
ential equations that a perturbation must satisfy. 


2 Linearized Perturbation Equations and 
Main Results 


To begin this section, we fi rst derive the system of PDFs 
that perturbations of (VP) must satisfy so that we may 


study their behavior. To determine the stability properties 
of these solutions we consider the initial value problem 
with time dependent perturbations of the form: 


fit, x,v) = p. -f (b{x) ) -f F{t, x, v) 
Eft, x) = (bfx) -I- j3{t, x). 


Using these quantities in (VP) yields equations for the 
perturbations, namely 


dtF -\- v(b'p! -\- vdxF — {(b' -I- l3)vp,' — {cf' -I- b3)dyF = 0 
4/' d- dxP = 1 — J fi{e)dv — J F{t, X, v)dv. 


Since (/r, 4>) is a time-independent solution to (VP), these 
equations simplify and we arrive at 


dtF -\- vd^F — (b'dyF = vjbg,' d- jbdyF' 
dxP = - [ Fdv. 


( 2 . 1 ) 


Noti ce that the system of partial differential equations 
(2.1) is nonlinear due to the appearance of the fbdyF term. 
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We focus on studying the linearized system, so we remove 
this term, which yields 


which is identical to the first equation of (2.3) under the 


dtF + vdxF — (jj'dyF = vPfj,'' 
dxP = - [ Fdv. 


( 2 . 2 ) 


Theorem 1 Let A S K 6 e given Then, equation (2.2) pos¬ 
sesses a solution of the form 


F{t,x,v) = e^*F{x,v) 
j3{t, x) = e^^'F{x) 

if and only if it possesses a solution of the form 

F{t, X, v) = e~^^F{x, —v) 
f3(t, x) = e~^^'F{x). 

Proof (Theorem^ The result of the theorem follows from 
the time-reversible Hamiltonian nature of (2.2) (see [12] 
for a more general result concerning the Vlasov equation). 
In this case, we merely utilize the specific forms of solu¬ 
tions to the linearized equation. Namely, using 

F{t,x,v) = e^*F(x,v) 
j3{t, x) = e^^F{x) 


within ( 2 . 2 ), we see that this is a solution if and only if 


\F{x, v) -\- vdxF{x, v) — (j)'dvF{x, v) = —if'vp!' 




(2.3) 


Next, we instead impose 

F{t, X, v) = e~^*F{x, —v) 
j3{t, x) = e~^^F{x) 


within ( 2 . 2 ) and find 


for unknown perturbations of both the distribution func¬ 
tion F{t, X, v) and the field I3{t, x). This linearized system 
of PDE then constitutes the equations under study. Thus, 
the main question of interest is whether solutions {F, j3) 
tend to zero over long times, thereby leading the solution 
{f,E) of (VP) to tend to the BGK waves, or whether so¬ 
lutions to may remain large over time, thereby caus¬ 
ing / and E to stay far from equilibrium. Our first result 
show s that exponentially growing and decaying modes of 
( 2 . 2 ) must occur in pairs of opposite sign. 


relabeling w = v. Changing variables using w = —v within 
the integral of th e second equation then yields the second 
equation of (2.3). Therefore, the two solutions satisfy the 
same system of equations, and the proof is complete. 

Therefore, one cannot expect asymptotic stability of 
solutions as all decaying modes must give rise to a grow¬ 
ing mode. We note, however, that Theorem ^ does not 
guarantee the existence of either solution, and this is the 
focus of the remainder of the paper. Our main result con¬ 
cerning the linear instability of periodic BGK waves (1.1) 
- ( 1 . 2 ), is stated precisely in the following theorem. 


Theorem 2 Let f and E he the periodic BGK solutions 
of {VP) _described by (1.1) with period P^. Assuming con¬ 
ditions ([^ - (Q, there is A > 0 and a solution of the form 

F{t, X, v) = e^*F{x, v) 

P{t, x) = e^^F{x) 

to the linearized perturbation equations [2^ where F{-,v) 
is a P^-neriodic function for every v G M, and F = —ib' 
wUh/GH\0,P,). 

Hence, there exists a growing mode for the linearized per¬ 
turbation equations with period Prf,, and from Theorem 
one can construct a corresponding decaying mode from 
this solution, as well. Many authors have investigated ques¬ 
tions in this vein under specific conditions on the func¬ 
tion q, for instance, see [1] and |13j for results concerning 
single-well and convex functions, respectively. In the case 
of Theorem the properties of q greatly depend on the 
local behavior of the distribution function p. and potential 
(j), and even for different single-well or convex profiles, as¬ 
sumption (Q may be satisfied or fail to hold. Regarding 
the proof, the main idea hinges on utilizing an equiva¬ 
lent formulation of the problem developed in m- Prior 
to stating the lemma containing this result, we first con¬ 
struct the characteristic curves for the particles in ( 2 . 2 ). 
Hence, define the curves X{s,x,v) and V{s,x,v), which 
we shall often abbreviate as X (s) and V (s) respectively, as 
solutions of the system of ordinary differential equations 


ds 

X(0, X, v) = X, 


= y{s,x,v), 


dV 

ds 


= -4>'{X{s,x,v)), 
V (0, X, v) = V. 


(2.4) 


—\F{x, —v) -\- vdxF{x, —v) -\- (j)'dyF{x, —v) = —if'vp' 
tj}"{x) = JF{x,—v) dv. 

Then, writing w = —v in the first equation and removing 
a negative sign in front of each term, we find 

\F{x, w) -\- wdxF{x, w) — 4> dyF[x, w) = —if'wp! 


From the definition of these curves, we can immediately 
construct an invariant of the physical system, namely the 
particle energy 

e=]^v‘^ + (j){x) = \v{sf + <j){X{s)) 
which is independent of the time variable s . Th is is ob¬ 


tained by multiplying the equation for ^ in (2.4) by V(s) 
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and expressing the remaining terms as a total derivative in 
s. Throughout, we will denote both the particle energy and 
the exponential function by e, leaving the reader to differ¬ 
entiate between them due to context. For instance, we will 
often write /r(e) in the future instead of ^ 
or fi {\V(sY + even though the three terms are 

equivalent. With this in hand, we may now summarize 
the previous work of m in the context of singly-periodic 
perturbations of BGK waves. 

Lemma 1 (in]) Let X > 0 be given. There exists a non¬ 
trivial solution to of the form 

F{t,x,v) = e^*'F{x,v) 
j3{t, x) = e^^'T{x) 


Figurej^contains a representation of the periodic wave 
(j) and displays the associated periods of the particles oc¬ 
curring within In view of Lemmawe can now see 
that the sign of fj! {e) may play a large role in the stability 
of solutions. For instance, if /i'(e) < 0 for e G [(/>_, oo) one 
could guarantee the last two expressions in £[■0] remain 
non-positive while choosing 0 to minimize the first inte¬ 
gral. However, as we will show, BGK waves cannot support 
such distribution functions, meaning /r^(e) cannot remain 
nonpositive for all energies e G [0_, oo). 


Theorem 3 Assume conditions |z|)-|zt|) hold and the BGK 
wave (/i, 0) satisfies (1.2), then 


q{x) 





dv > 0 


if the functional £[0] defined by 


r 


£[ 0 ] := 


-k2 


|0'(x)|^ — y fi'{e) dv\if{x)\'^ dx 

4+ Pf \Jo \/2(e - (j){y)) 

h- ^ Pt \Ja \j2{e - 


dy ) de 


dy 


for some x G [0, P^]. In particular, p cannot satisfy p'{e) < 
0 for all e G [0_, oo). 

We remark that the second statement of the the theorem 
regarding the non-monotonicity of /i(e) has been shown 
before within and [3S]. Hence, the signs of the com¬ 
plicated terms in the representation of Lemma depend 
very delicately on the behavior of p'{e). Regardless, we 
dc. can eliminate this difficulty by utilizing the BGK waves 
and their symmetry to prove Theorem 


satisfies £[0] < 0 for some 0 G H^{0,P^). Here, a G 
[0, lP<t>] is the unique point in the interval such that 4>{a) = 
e, for a given e G [0_, 0+], while 


Pf = 


1 


--dy, 


0 \/2(e- 0(y)) 

X 

Ja \/2(e- 0(j/)) 

Additionally, J- and W are given by 

P{x,v) = p'{e) i^j Ae'''®0(Al(s))(is — 0(a;)^ (2.5) 


and F = —if' , where where X{s) is the spatial solution of 
the characteristic equations [2^. 

We refer the reader to [T^ for the details of the proof, but 
we can also briefly sketch the main ideas. More specifically, 
utilizing the form of a growing mode solution (F, /?) from 
Theoremj^within the linearized perturbation e quati ons, it 
is a straightforward computation to show that ( | 2 . 2 [ ) has a 
solution if and only if Axif = 0 for an associated family of 
Hermitian dispersion operators A\. With this, m proves 
the continuity of A\ with respect to A > 0 and positivity of 
Ax as A —> -koo, while the condition of Lemma namely 
£[0] < 0, guarantees that Ax is negative as A —>■ O’*". 
Hence, an application of the Intermediate Value Theorem 
then implies the existence of a value of A > 0 for which 
Axip = 0. To keep technical details to a minimum, we will 
omit a complete description of the dispersion operator Ax 
and focus instead on proving the condition £[0] < 0 . 


Proof (Theorem^ Since the latter terms in the represen¬ 
tation of £[0] are complicated, we will construct a specific 
choice of 0 G H^(0,P^) for which the first piece is neg¬ 
ative and use symmetry to eliminate the contribution of 
the last two expressions. To begin, recall (1.2): 


4>"{x) = 1 — J p{e)dv. 


We take an cc-derivative to find 


0"'(x) 


p'{e) dv 


(j)'{x) 


or 

(j)'"{x) + q{x)(j)'{x) = 0 . 

Hence, the function u = (f>' satisfies the Sturm-Liouville 
problem 

u"{x) -k {q{x) -k X)u{x) = 0 

u(0) = u{P^) = 0 

for u G H^{0,Prj,) with corresponding eigenvalue A = 0. 
However, 0 attains a maximum at and possesses no 
other critical points on (OjP^). Thus, (j)'{^P^) = 0 and 
because 0 ' has exactly one root in (O,P 0 ), it must be the 
second eigenfunction [ 23 . Therefore, there exists a first 
eigenfunction uq G H and eigenvalue Aq < 0 such that 

Mo(a:)-k (g(x)-k Ao)Mo(a:^) = 0 (2.6) 

and Uo(0) = Uo{P,p) = 0. Using the symmetry of q, namely 
qiP^, — x) = q{x) for every x G [0,P,j,], we see that the 
function uo(x) := uo(P 0 — x) is also a solution of the 
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Fig. 1 : A representative graph of the potential 0(a;) for the BGK waves of ( |l.l[ )-( |l.2| |. The particles are divided into free 
particles (with period Pf), whose energy lies completely above the dashed line (e > 0+) and hence above the graph of (j), and 
trapped particles (with period Pt), whose energy lies on or below the dashed line {(j>- < e < <j>+) and intersects the graph of 
(j). For the latter, the symmetry of cji implies that the dotted line denoting the particle energy must intersect the graph of ij) in 
exactly two points within the interval [0, P,^]. One of these points (denoted by a) must be less than \Ptj, and the other, at the 
reflected point P^ — a, must be greater than |P,^. 


Sturm-Liouville problem with the same eigenvalue. Due to 
uniqueness of solutions, it follows that there is C S K such 
that uo{x) = Cuo{P^ — x). Evaluating this at the point 
X = ^P4>, we see that either C = 1 or uo{^P^) = 0. As uq 
is the first eigenfunction, it has no roots in {0,P^). Thus, 
(7 = 1 and uo(x) = uo{P^ — x). Since mq G H‘^(fi,P^) 
we may take a derivative of ( |2.6[ ) and find that Uq sat¬ 
isfies the symmetry condition Uq{P^ — x) = —Uq{x) for 
every x € [0, P,^]. With these functions in place, define 
ipix) = uo{x)uq{x) for X G Using the aforemen¬ 

tioned properties of uq and Uq, we see that Tp satisfies 

V^(O) = V'(P0) = 0 

with ip{P^ — x) = —ip{x) for all x G [0, P^]. In addition, ip 
can be extended to the whole space in a smooth manner 
by imposing P^-periodicity since 

Pj'iP^) = |u'(P^)|2+uo(P0K(P0) 

= 

= l^o( 0 )r 

= K(0)p + uo(0K(0) 

= ^'{0)- 

As Mo G we see that ip G H^{0,P^) and is 

smooth at the endpoints. Furthermore, this function sat¬ 
isfies 


\P:f-q\Pj\^ = \{u'of+Uou''\^-q\uou'o\^ 

= |K)^ - (<? + Ao)(mo)^|^ - q\uouX 

= kol‘‘-2(g-bAo)|MoMoP 
+(9 + Ao)^|uo|‘* — q\uQUQ\^ 

= ImqI"* — 2Ao(moMo)^ + (9 + Ao)^|mo|^ 
-3g|uoMoP 


Finally, we compute each portion of the representation 
given in Lemma For the first integral, we have 
rPd. r 


- q{x)\i:{x)p 


dx 


i-Pa V 


uX - 2Ao|moUo 7 -b (9 + Ao) |mo| - 39 |moMo|" 


dx 


For the first term of the integral, we integrate by parts 
and use (2.61 to yield 

f XoixX dx = f (^uo{x)] {uq{xX dx 


= -3 


uo(x)|Mo(a;)pMo(a:) dx 


= 3 


{q{x) + Ao)|uo(a;)Mo(a;)|^ dx 


Notice here that we have utilized the conditions Mo(0) = 
Uq {PcfP) = 0 to eliminate boundary terms arising from in¬ 
tegration by parts. Hence, we find 


XXX - <i{xM{xX 


dx 


AoI^oMqP -b (9 -b Ao)^|uo|^ 


dx. 


The first term o f th is expression can be simplified further 
by multiplying (2.6) by Mg and integrating so that 

i-P^ f-P^, 

/ {q{x) -b Ao)|Mo(a;)7 dx = - UQ(x)(uo(x)f dx 

Jo Jo 

= ( [(wo(a:))^] dx 


= 3 


|Mo(a:)Mo(a;)|^ dx. 
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Using this along with the even symmetry about of 
|mo|"^ and q, the first piece of the representation satisfies 

r 1 




q{x)\ii){x)\^ 


dx 


Jo 


fPd, r 


2 ^ 0(9 + Ao)|uo|"^ + {q + Ao)^|uo|^ 


(g + Ao) ( g + -Ao ) |uo| dx 


dx 


Now, since ^|J is odd about |_P 0 , the remaining terms in 
the representation of LemmalRwill vanish. More precisely, 
for every e G [^+, 00 ) we have 



\/ 2 (e- (^(y)) 


dy = 0 


(2.7) 


and for every e G [(()_, ^+] we have 



tpjy) 

y/2{e-(j){y)) 


dy = 0. 


( 2 . 8 ) 


To prove (2.7), we first decompose the integral into the 
reflected components of tp, 


r , 'I'j 

Jo \/ 2 (e - (piy)) 


r ‘ -PM=dy 

Jo \/ 2 (e - </>(?/)) 

+ /'■* -^tM}=dy 

JhPd> y/2(e - ((-(y)) 

I+ 11. 


Now, computing term II, we change variables z = — y 

and use the symmetry of ip and (p from ( p^ to demonstrate 
the cancellation of these terms, so that 

„ = r 


hp* \/ 2 (e - ())(y)) 

^(z) 


lo y/2{e- (p{P^ - z)) 


2 P^ 


= -I 


dz 


i;{z) 


lo ^2{e-(p{z)) 


’'2^4' 


ip{z) 


lo dy2{e-(p{z)) 


dz 


dz 


Hence, adding thes e terms yields (2.7). Turning to the 
justification of ( 2 . 8 ), we again divide the interval of in¬ 
tegration about the midpoint Using the symmetry 

and periodicity of (p and u, we perform the same change 


of variables as above to find 

4iy) 


dy 


hPd. v/2(e- (()(y)) 

Pd{P^ - y) 


'hPd, \/ 2 (e - </>(y)) 


dy 


"pPd> 


4{z) 


dj2{e-(p{P^ - z)) 


dz 


'■\P4> 


-Ipiz) 


y/2(e-(p(z)) 


dz. 


Hence, we conclude (2.8) by adding these terms. 

Finally, from these computations and in view of as¬ 
sumption (T^, we find 


Clip] = 


\4'{x)\ - / y'(e) dv\pj{x)p 


dx 


+ M'(e)yw 


V'(y) 


0 + Pf \Jo ^J2{e-(p{y)) 

4{y) 


dy ) de 


r<i>+ I 

+2 / y'(e)- 

Jd- 


cP(h — Oi 


\/ 2 (e - (piy)) 


dy ) de 


< 0 . 


~pPd> 


{q{x) + Ao) ( q{x) + ^Ao ) |uo(a:)|'‘ dx 


which proves the main result. 

Finally, we conclude with the proof of Theorem 
Proof (Theorem^ Let /r G and (p G C'^ satisfying 

(p" {x) = 1 - J y. + (p{x)^ dv (2.9) 

and assumptions (|^ - ( p^ be given. We will take advantage 
of the regularity ot thes e fu nctions. Assume q{x) < 0 for 
all X G [ 0 , iP^] . Using (2.9), we see that the regularity o:’ 
the right side implies that (p G and 

(P'"{x) = -q{x) (p'{x). 


Because g assumes only negative values, we see that (p"’{x) 
and (p'{x) must have the same sign for all x G [O, jP^]. 
From (iv), (p is strictly decreasing on the interval [0, ^P^] 
with a maximum and minimum at ^(0) and (p (jP^), re¬ 
spectively. Thus, </>'(0) = (p'{^Prj,) = 0 and (p'{x) < 0 
on the interval (O, ^P^)- Therefore, (p'" < 0 on (O, |P(/,) 
and thus (p” is decreasing on (O, \P^)- However, from the 
above conditions on (p', we see that cp' must transition from 
decreasing to increasing on some subinterval of (O, ^P^), 
and this implies that cp” must transition from negative 
values to positive values on the same subinterval. Clearly 
such a transition cannot occur if (p"{x) is decreasing for ev¬ 
ery X G ( 0 , ^P^) and we arrive at a contradiction. There¬ 
fore, there exists x G [O, ^P^] such that 

q{x) = J d' dv > 0. 
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Additionally, this condition cannot be satisfied for any /i 
with fi'{e) < 0 for all e G oo) and this completes the 
proof. 

3 Conclusion 

In summary, specific conditions were determined that guar¬ 
antee the linear instability of periodic BGK waves for the 
one-dimensional Vlasov-Poisson system. This result then 
limits the variety of BGK waves which may be stable un¬ 
der singly-periodic perturbations. In particular, our main 
results demonstrate that for a BGK wave to be stable, 
specific values of a quantity involving the derivative of 
the particle distribution function cannot be essentially lo¬ 
calized within a given interval. While this sheds some light 
on the nature of BGK waves, it remains unknown as to 
whether they are stable or unstable with respect to general 
perturbations which share their period. Hence, it seems 
that the stability properties of these solutions greatly de¬ 
pend on the local behavior of the distribution function and 
potential. Additionally, it was shown that modes of the 
linearized system occur in pairs of opposite sign, thereby 
displaying that for every perturbation that decays expo¬ 
nentially in time there must be a corresponding pertur¬ 
bation that grows exponentially in time. Finally, it was 
shown that the particle distribution function cannot be 
monotonically decreasing for all energy values, which also 
limits the stability properties of these waves. 
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